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Abstract 



Cf^ ' A quantum master equation is obtained for identical fermions by including a relaxation term in 

■^ ■ addition to the mean-field Hamiltonian. [Huang C F and Huang K N 2004 Chinese J. Phys. 42 

in ■ 

O ; 221; Gebauer R and Car R 2004 Phys. Rev. B 70 125324] It is proven in this paper that both 



the positivity and Pauli's exclusion principle are preserved under this equation when there exists 
an upper bound for the transition rate. Such an equation can be generalized to model BCS-type 
Mh, quasiparticles, and is reduced to a Markoff master equation of Lindblad form in the low-density 



Ch ' limit with respect to particles or holes. 



PACS numbers: 



I. INTRODUCTION 



The quantum kinetic approaches have attracted much attention because of the develop- 
ment of the nanoscale devices, for which the classical kinetic theory may become invalid. 
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3een introduced to develope the quantum 
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17j In these equations, there are 
relaxation terms responsible for the irresversibility. The Markoff master equation p, [ij] is 
used to describe quantum relaxation processes, and the master equation of Lindbald form 
Iql is derived based on mathematical assumptions. After ignoring the pure-dephasing 



m 



term 



14l |. the Markoff master equation is just a particular master equation of Lindblad form. 



The well-known Pauli master equation, in fact, can be deduced from the Markoff master 
equation in the incoherent limit. IJ] 

It is important to assume the linearity to derive the master equation of Lindblad form 
in quantum mechanics, p^] On the other hand, nonlinear relaxation terms have been intro- 
duced semiclassically based on Pauli's exclusion principle for systems composed of identical 
fermions. [l8|, [l9(] A nonlinear quantum master equation has been introduced in Refs. 8] 
and [lOy to unify the quantum and semiclassical approaches. The relaxation term of such 



an equation, in fact, can be constructed by considering two antihermitian terms. [Sj] One is 
responsible for the loss of particles while the other is for the loss of holes, which is equiv- 
alent to the gain of particles. Here holes are vacencies of any orbitals. (8|, [l9(] Therefore, 
the relaxation term is symmetric with respect to particles and holes. After incorporating a 
pairing tensor, the density matrices for particles and holes can be used to constructed those 



, which include 



23 



2J] have 



for quasiparticles in Bardeen-Cooper-Schrieffer (BCS) pairing models 

many-body effects beyond one-body approximation. Different approaches [22|, 

been developed to study many-body correlations, and BCS pairing models are powerful to 

understand superconductivity [25], superfiuid |2a], and meson-nucleon couplings 20l ]. 



It wi 
and 



1 be proven in this paper that the nonlinear master equation introduced by Refs. 
lO] preserves both the positivity and Pauli's exclusion principle when there exists 



an upper bound for the transition rate. Therefore, such an equation is suitable to model 
fermions. For convenience, first we discuss different types of master equations in section II, 
and the proof is in section III. As discussed in section IV, the nonlinear master equation 
for fermions can be reduced to the Markoff master equation in the low-density limit with 



respect to particles or holes. It is discussed how to prove the conservation of the trace of the 
density matrix. An extension of the master equation to the relativisitic Hartree-Bogoliubov 
model 



20 



2l| . which is a BCS-pairing model, is obtained by introducing a constraint on the 



relaxation term. It is also mentioned in section IV how to consider the re-pairing between 
BCS-type quasiparticles and quasiholes to incorporate multiple order parameters [26|, |28 |. 
We note that multiple order parameters have been introduced not only to understand heavy 
fermion metals [28j, but also to unify BCS theory and ferromagnetic/antiferromagnetic 



theory to probe high-temperature superconductors. 29|| Conclusions are made in section V. 



II. DIFFERENT TYPES OF MASTER EQUATIONS 



cesses 



lerent types ol 



Different types of master eq uations have been introduced to describe irreversible pro- 
lo|, y, y, y, y, y, y, y, y, y, y For & quantum system with 



the density matrix p(t), the master equation of the following form 

^^p{t) = i[p{t),H{t)] + R{p{t)) 



(1) 



has been discussed in the literature. Here R{p{t)) is the relaxation term, [v4, i?] denotes 
the commutator of any two operators A and B, and H{t) is the Hamiltonian generating a 
unitary operator U(t) by 



t^Uit) = H{t)U{t) 



(2) 



with f/(to) = I) the identity operator, at the initial time tg- In this paper, we take the 
reduced Planck constant /i = 1, denote Vt^ as the adjoint of an operator fi, and assume that 



all the kets are in a separable Hilbert space. [30|, l3l| Let |||«)|| = (a|a)^/^ for any |a), and 
we require that limt^^t^\\[U{t2) -U{ti)\\^)\\ = limt,^^t^\\[W{t2) -U^{ti)]\f3)\\ = for all |/3). 
We take the Hamiltonian H{t) = Hq + V{t) with two self-adjoint operators Hq and V{t) 
as the (time-independent) unperturbed and (time-dependent) perturbed parts. In addition, 
assume that Hq can be diagonalized by orthornormal complete basis such that each eigenket 
satisfying 



Ho\n) = En\n) 



(3) 



and \\\n)\\ = 1 can be parametrized by a positive integer n. Here each En is an eigenvalue. 



For quantum optics, 



14| the relaxation term is composed of the pure-dephasing term 



RMpipit)) = - ^ -fnn'\n){n\p{t)\n'){n'\ 

(n'n) 



(4) 



and the transition term 

RMt{p{t)) = "2 5Z ^nn'{p{t), W){n'\} + ^ Wnn'\n){n'\p{t)\n'){n\. (5) 

(n'n) (n'n) 

Here each orbital in Eqs. (4) and (5) is a (normalized) eigenket of Hq introduced by Eq. 
(3), Wnn' is the nonnegative coefficient for the transition from ket \n') to another ket \n) if 
n 7^ n', '^nn' is the nonnegative coefficient for the pure-dephasing rate of {n\p{t)\n') when 
n 7^ n', and we denote {A, B} = AB + BA for any two operators A and B. For convenience, 



we take w. 



nn inn 



In 



for all integers n in this paper, so we do not need to set the condition 



n ^ n' in the summations over {n'n) in Eqs. (4) and (5). If there is no pure-dephasing term, 
we can reduce the Markoff master equation as 



^^p{t)=i[p{t),H{t)]+RMt{p{t)) 



(6) 



i[p{t),H{t)] - - ^Wnn'{p{t),\n'){n\]+ ^Wnn'\n){n\p{t)\n')(n\. 



(n'n) 



(n'n) 



On the other hand, the famous relaxation term of Lindblad form 15|, [iGj 

RL{p{t)) = -\Y.{p{t)yM} + Y.^lp{m 
I I 

has been derived mathematically, where {V;} is a set of operators. By setting 

/ = (n'n) and V/ = wj^,\n'){n\, 



(7) 



(8) 



we can re-obtain RMt{p{t)) from Rii^pit)). Therefore, Eq. (6) is a particular Markoff master 
equation of Lindblad form. 

The relaxation term Rutipit)) in Eq- (6), in fact, is composed of the loss factor 



LM{p{t)) = -:^ XI Wnn'{p{t), \n'){n\} 



(9) 



(n'n) 



and the gain factor 



GM{p{t)) = y^ Wnn' \n) {n'\p{t) \n') (n| . 

(n'n) 



(10) 



With some calculations, we can see that the loss (gain) factor is to decrease (increase) the 
number of particles. Define the function f{n,t) = {n\p{t)\n), the loss and gain rates are 
Yl,n''^ri'nf{n.,t) and 'Ylm''^'nn'f{n',t) for cach orbital |n), respectively, [ij] In addition, from 
Eq. (6) we have 



d_ 
di 



{n\p{t)\n') = i{En' - En){n\p{t)\n) - - ^(wm„' + Wmn){n\p{t)\n') 



(11) 



for all phases {n\p{t)\n') with n 7^ n' if the perturbed potential V{t) = 0. Thus the Markoff 
master equation can incorporate decoherent effects. [l3|, [ij] The decoherence comes from 
the loss factor LM^pit))- From Cauchy-Schwartz inequality, we have 



\{n\p{t)\n')\< {n\p{t)\nY'^n'\p{t)\n') 



.'\l/2 



(12) 



If {n\p{t)\n) —>■ because of the loss effects on orbital \n), the phase {n\p{t)\n') should 
also approach zero for any other n' based on the above equation. Since Cauchy-Schwartz 
inequality is valid for all positive matrices, it is natural that the loss term results in the 
decoherence to preserve the positivity of p(t). When H{t) = Hq, actually the Pauli master 



equation 



i3i,y,ii8| 



d_ 
di 



f{n, t) = -^ Wn'nfin, t) + ^ Wnn'fin', t) 



(13) 



can be obtained from Eq. (6). 

For a system composed of noninteracting identical fermions, usually the matrix p{t) is 
taken as the one-particle density matrix with the trace trp{t) = X]n('^lp(^)l^) ~ X]n/(^'^) 



equals the number of particles. 



12j In such a case, the function f{n,t) = {n\p{t)\n) is 



interpreted as the occupation number in orbital \n) at time t. It is known that Eq. (13) 
should be modified for fermions to follow Pauli's exclusion principle. 18| By extending the 
coefficient Wnn' as a time-dependent function 



Wnn'(t) = LJnn'[l " f{n,t)], 



the master equation of the following form 



d 



aisi. 
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18. 


19 



-f{n, t) = -Y, ^n'«[l - /(^', i)]/(^, t) + Yl ^««' [1 - •/'(^' ^)]-/'(^'' ^) 



at 



(14) 



(15) 



can be obtained for fermions from Eq. (13). Here Unn' is the nonnegative coefficient for the 
transition from n' to n when n 7^ n', and we set Unn = for all n. On the other hand, the 



nonlinear quantum master equation 



^p(t) = iW),H{t)] -]^Y. ^-'(1 - {n\p{t)\n)){p{t), \n'){n'\} (16) 



{n'n) 



+ \ Y. ^nn'{n'\p{t)\n'){I - Pit), \n){n\} 



2 

(n'n) 

is introduced recently for systems composed of noninteracting identical fermions. [8| The 
operator 

p(^)(t)=/-p(t), (17) 



which appears in the last term of Eq. (16), can represent holes actually. [Sl, loj In Ref. 
Eq. (16) is obtained from the following equation 

^p(t) = t[pit), Hit)] + {pit), A,it)} -{I- Pit), A^it)} (18) 

by considering the conservation of the number of particles in each transition. If we define 

A,in) = IY1 ^""'(1 - {n\n\n))\n'){n'\ (19) 

{n'n) 

v42(f2) = 2 5Z ^nn'{n'mn')\n){n\, 

(n'n) 

for any operator Q, Eq. (16) can be obtained from Eq. (18) by setting Apit) = — Ai(p(t)) 
and Apit) = — A2(p(t)). We can see that Eq. (16) is equivalent to the one-body master 



equation introduced by Ref. [lOj after expanding p(t) with respect to the eigenorbitals of 
Hq. Such an equation is denoted as Kohn-Sham master equation in Ref. p] because it can 
be used to extend Kohn-Sham equations. When H(t) = Ho, Eq. (15) can be derived from 

fin 

Eq. (16). [8|, 110[ On the other hand, we can reduce Eq. (16) to Eq. (6) in the low-density 
limit by setting u;„„/ = Wnn'- [8] Therefore, Eq. (16) can be used to unify Eqs. (6) and (15). 
The second term at the right hand side of Eq. (16) 

Lfipit)) = -^ E ^-'(1 - {n\pit)\n)){pit), \n'){n'\} = -{p(t), Ai(p(t))}, (20) 

(n'n) 

in fact, can be obtained from the loss factor L^fipit)) of Eq. (6) by Eq. (14). With some 
calculations, it is easy to see that Lfipit)) induces the loss of particles just as LA/(p(t)) does. 
Hence Lfipit)) serves as the loss factor in Eq. (16). For each orbital |n), the loss rate due 



to Lf{p{t)) equals Yln'^n'nl^ — f{n',t)]f{n,t). In addition, the factor Lf(t) also results in 
the decoherence, which is important to preserve the positivity of p{t) as mentioned above. 
On the other hand, the third term of Eq. (16) 

GfiPit)) = ^ E ^nn'{n'\pit)W){I - Pit), \n){n\} = {I - p(t), ^(p(t))} (21) 

(n'n) 

induces the gain of particles and serves as the gain factor. For each orbital \n), the gain 
rate due to Gf{p{t)) equals J^n'^nn'l'^ — f{n,t)]f{n',t). Such a factor, however, cannot be 
obtained from the gain factor GM{p(t)) of Eq. (6) by Eq. (14). With some calculations, we 
note that Gf{p{t)) results in the decoherence while GM{p(t)) does not. 

The gain of particles (holes), in fact, is just the loss of holes (particles) in each specific 
orbital. Sj Hence it is natural to relate the the gain factor Gf{p{t)), which can be taken as 
the loss factor for holes, to Lj{p{t)). We note 

LMt)) = -\Y. ^nn'{n\p^\t)\n){I - p^\t\ \n'){n'\], (22) 

(n'n) 

from which we can obtain the gain factor Gf{p(t)) by Unn' — ^ ^n'n , changing the sign, and 
relacing p^^^'it) by pit). The process that particles jump from n to n' should correspond 
to the transition of holes from n' to n, and thus it is reasonable to replace Unn' by Un'n 
to construct Gf{p(t)) from Lf{p(t)). The change of the sign is natural because the loss 
factor is to decrease the occupation number while the gain factor is to increase it. We shall 
replace p®(t) by p(t) since Lf{p{t)) and Gf{p{t)) correspond to the gain of holes and that 
of particles, respectively. 

As mentioned above, the decoherence due to the loss factor Lf{p(t)) (for particles) is 
important to preserve the positivity of p{t). Since Gf{p{t)) can be taken as the loss factor 
for holes, it is natural for Gf{p{t)) to result in the decoherence to preserve the positivity of 
p^p\t). The positivity of p^\t), in fact, is equivalent to that p{t) follows Pauli's exclusion 
principle because {a\p{t)\a) < 1 iff 

{a\p^^\t)\a) = 1 - {a\p{t)\a) > (23) 

for any normalized \a). Therefore, the decoherence is important not only to the positivity, 
but also to Pauli's exclusion principle. 



III. THE POSITIVITY AND PAULI'S EXCLUSION PRINCIPLE 

To show that Eq. (16) is suitable for fermions, it will be proven in this section that both 
the positivity and Pauli's exclusion principle are preserved when there exists a positive real 
number M such that both 

sup ^ -Wnn' and sup ^ -Wn'n < M. (24) 

n' n' 

For a specific \n), the loss and gain rates J2n''^n'n[^ — f{^'{t)]f{^,t) and J2n''^nn'[^ — 
f{n, t)]f{n', t) should be both smaller than 2M. So such a number provides an upper bound 
for the transition rate. For convenience, let to be the initial time for the time evolution 
of pit). Just as mentioned in the last section, the unitary operator U{t) generated by the 
Hamiltonian H{t) equals the identity operator J at t = to- We require p(to) is a positive 
self-adjoint operator following Pauli's exclusion principle such that 

< (a|p(to)|a) < 1 for all normalized |a). (25) 

Because Eq. (16) yields an initial-value problem, we just need to prove the existence and 
uiqueness of the solution when t is in a time interval [to, tj] for some final time tj satisfying 

to<tf<t, + ^. (26) 

In the following, the notation \\Q\\ = sup|||„\||^;^ ||(5|tt)|| as the natural norm of any bounded 



(linear) operator Q |30|, l3l|- It is known that ||(5^|| = \\Q\\ under such a norm. 3l| In 
addition, we denote |||Q(t)||| = sup^gj^^^^^ i ||Q(t)|| for any time-dependent operator Q(t) if 

suPte[to,t^] 1 1 Q{t) 1 1 is finite. 

It is convenient to introduce the following Banach spaces 30|, l3l|, l3^ Si and S[, a subset 
^2 of Si, and a vector space Sq: 

Definition 3.1 Let S'l be the Banach space composed of all the mappings Q'{t) from the 
time interval [to,t/] to bounded (linear) operators such that 

lim ||fi'(t2)-fi'(ti)|| = 0, (27) 

and denote the Banach space Si = {Q{t)\Q{t) = U{t)Q'{t)U'^{t) for some Q'{t) G S[}. 
We take \\\ ■ \\\ as the norm on Si and S'l- The set S2 = {^{t) G Si\fl{t) = f2'l'(t),0 < 



(Q;|n(t)|Q;) < 1 for all normalized \a)} is a complete subset of Si. The vector space So 
is composed of all the mappings Q{t) from [to,t/] to bounded operators such that 

|||Q(t)||| < ooand lim ||[Q(t2) - Q(ti)]|a)|| = for all |a). (28) 



The unitary operator U(t) and its adjoint U^(t) are in Sq. A density matrix p(t) (in 
the Schrodinger picture) is a postive one following Pauli's exclusion principle if p{t) G S2- 
For each operator Q(t) G Si in the Schrodinger picture, its corresponding operator in the 
Heisenberg picture is Q'(t) = W(t)Q{t)U{t) G S[. The adjoint of any operator in S[ is 
also in S'l. Actually Si and S'l are subspaces of So- For convenience, any operator in 
S[ is denoted by the superscript "/" in this paper. As discussed in Appendix A, we can 
denote the operator O = j].^dt'Q{t') for any Q{t) E Sq iS O is the uique one following 
{a\0\f3) = j'l'^ dt' {a\Q{t')\l3) for all \a) and \I3). Here ti and ts follow to < ti < ta < tf. In 
addition, we have 

\\0\\ = \\ I' dt'Q{t')\\<\t2-ti\x\\\Q{t)\\\. (29) 

Jtx 

For any time-dependent operator Q{t) G iSq, we denote 

|-fi'(t) = Q{t) for some 9.'{t) G S[ (30) 

on an interval [^1,^2] ^ [^o,t/] iff 

n\t) = n\ti) + f dt'Q{t') (31) 

at any t G [ti,t2]. It should be noted that such a definition can be invalid for §[^{t) when 
Q{t) G Si. Some properties about the integral and drivative of operators are discussed in 
Appendix A. 

The following lemma is important to the proof: 
Lemma 3.2 Consider the equation 

-p{t) = i[p{t),H{t)] - {p{t),A{t)} + B{t) as t G [to,t/], (32) 

and assume that the given initial matrix p{to) satisfies Eq. (25). Here A{t) and B{t) 
are two positive self -adjoint operators in Sq, and we assume that |||v4.(t)||| < 2M. Then 



there exists a unique self -adjoint matrix p{t) G Si to follow the above equation, and the 
positivity is preserved. 

proof To avoid unbounded problems due to the Hamiltonian H{t), we can rewrite Eq. 
(32) as 

^/{t) = -{p'it),Aum + Buit) (33) 

by introducing p'{t) = U\t)p{t)U{t) as the density matrix in the Heisenberg picture. Here 
Au{t) = W{t)A{t)U{t) and Bu{t) = W{t)B{t)U{t). We have p'(to) = p{to) at the initial 
time to because f/(to) = I- The matrix p'{t) G S'l iff p{t) G Si, and we can prove that both 
Au{t) and Bi/it) G Sq from corollary Al. Let JF be the mapping from S[ to S[ such that 
A'(t) = J^(fi'(t)) iff 

A'(t) = p(to) + tdt\Bu{t') - {n\t'),Au{t')}). (34) 

Jto 

The mapping JF is well-defined and bounded because we can use corollary Al to prove that 
the integrand Bu(t) — {Q'{t),Au(t)} G Sq. In addition, we can see the equivalence between 
Eq. (33) and p'{t) = J^{p'{t)). Since |||^c/(t)||| = 111^(^)111 < 2M, JF is a contraction on 
<S( and there exists a unique solution p'{t) G S[ to Eq. (33) from the fixed point theory. 



32l | We can take Eqs. (32) and (33) as the same equation in different pictures, so p{t) = 

U{t)p'it)W(t) is the unique solution to Eq. (32). 

Define a mapping K'{t; t') from t' G [t^, tf\ to S[ by 

d 

TrK'it; t') = -Au{t)K'{t; t') att>t' (35) 

and K'{t; t') = I a.t t < t'. In addition, let K{t; t') be a two-parameter operator such that 
K{t;t') = U{t)K'{t;t')U\t') when t > t' and K{t;t') = I when t < t'. It is shown in 
Appendix B that K'{t; t') can also be taken as a mapping from t G [to, tj] to S[, and we have 
§lK'\t; t') = -K'^{t; t')Auit) at t > t' from corollary A5. On the other hand, f if (t; t') may 
suffer unbounded problems if H(t) is not bounded although it is convenient to introduce the 
following equation 

i— K(t; t') = [H{t) - iAit)]K{t; f) at t > t'. (36) 

To prove the self-adjointness and positivity, we note that the solutions to Eqs. (33) and (32) 
are 

p'{t) = K\t- to)pito)K'^it; to) + / dt'K'it; t')Bu{t')K'^{t- 1') (37) 

J to 

10 



pit) = Kit- to)pito)K^it; to) + / dt'Kit; t')Bit')K\t- 1'), (38) 

J to 

respectively. It is discussed in Appendix B how to prove that both Eqs. (37) and (38) 
are well-defined. The above two equations can be related by p{t) = U {i) p' it)U'^ {i) , and we 
just need to check p'(^o) and ^p'{t) to see that they provide the solutions. We can use 
corollaries A2 and A4 to perform the time-derivative on p'{t). It is easy to see that p{t) and 
p'{t) constructed by the above two equations are self-adjoint. To prove that p{t) is positive 
such that ip.\p{t)\a) > for all \a), we note that 

{a\p{t)\a) = {a{t,toMto)\a{t,to)) + f dt'{a{t,t')\B{t')\a{t,t')) > (39) 

J to 

from Eq. (38) if we define \a{ti,t2)) = K^{ti;t2)\a). QED 



It is mentioned in the last section that Eq. (16) can be obtained from Eq. (18) by 
setting Apit) = —Ai{p{t)) and Ap(t) = — A2(p(t)). Thus it is convenient to discuss Eq. (18) 
before completing the proof for Eq. (16). Based on lemma 3.2, it is proven in the following 
proposition that Eq. (18) preserves both the positivity and Pauli's exclusion principle. 
Proposition 3.3 Let Ap{t) and Apit) be two time- dependent self-adjoint operators in 
Sq. Assume that both —Apit) and —Apit) are positive and \\\Apit) + Ap(t)||| < 2M. Then 
there exists a unique matrix pit) G ^2 satisfying Eq. (18) if Eq. (25) is valid. 
proof We can rewrite Eq. (18) as 

^p(t) = t[pit), Hit)] + {pit), Apit) + Apit)} - 2Apit). (40) 

By setting Ait) = —Apit) — Apit) and Bit) = —2 Apit), the above equation is of the same 
form as that of Eq. (32) and the condition |||.4.(t)||| < 2M is satisfied. Hence there exists a 
unique self-adjoint solution p(t) G Si following 

{a\pit)\a)>0 (41) 

for all normalized \a) from lemma 3.2. On the other hand, we can also rewrite Eq. (18) 
with respect to p®(t) = I — pit) as 

|p® (t) = ^[p^^\t), Hit)] + {p(^)(t), Apit) + Apit)} - 2Apit) (42) 



11 



The above equation is also of the form of Eq. (32), and Eq. (25) indicates that p*-^-*(t) is 
positive. Then we have 

1 - {a\p{t)\a) = {a\p^P\t)\a) > (43) 

for all normalized \a). We can complete the proof from Eqs. (41) and (43). QED 

Just as in lemma 3.2, we can rewrite Eq. (40) with respect to p'{t) = U'^{t)p{t)U{t) to 
avoid the unobunded problem due to H{t). In the Heisenburg picture, the density matrices 
for particles and holes are p'{t) and I — p'{t) = W (t) p^'P\t)U (t) , respectively. We can see the 
meaning of %p^'^\ which appears in Eq. (42), by considering the time- derivative on I — p'{t). 

Now we return to discuss Eq. (16). As shown in Appendix C, both Ai{VL{t)) and A2{^{t)) 
are positive self-adjoint operators in Sq for any Vt{t) G ^2. In addition, direct calculations 
yield that 

\\\A,{n^{t)) - A,(n,{t))\\\ < Mm,{t) - n,(t)\\\ (44) 

for any two operators ^i{t) and VL2{t) in ^2, where J = 1 and 2. The solution to Eq. (16), 
in fact, satisfies p{t) = F{p{t)) if we define the mapping F from ^2 to ^2 as: 
Definition 3.4 For any n(t) E S2, let 

F{n{t)) = JCnit; to)p(to)4(t; to) + 2 / dt'ICn{t; t')^(^(t'))4(^; ^')- (45) 

Jto 

Here lC^{t;t') follows 

2^/Cn(t; t') = [H{t) - iAi{n{t)) - iA2{n{t))\lCn{t- 1') at t > t' (46) 

and )Cn(t; t') = I at t < t'. Then 

^A(t) = z[A{t), H{t)] - {A(t), A,{Q{t))} + {I- A(t), A^m))} (47) 

= t[A{t), H{t)] - {A(t), Ai(fi(t)) + A2{n{t))} + 2A2{n{t)) at t > to 
and A(to) = p(to) iff A(t) = F(fi(t)). 

We can see the equivalence between Eqs. (45) and (47) by comparing them to Eqs. (38) 
and (32). While Eq. (46) may suffer the unbounded problem due to the Hamiltonian H{t), 
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we note that ]Cn(t;t') is just the solution to Eq. (36), which is equivalent to Eq. (35), 
if we take A{t) = Ai{n{t)) + v42(fi(t)) and B{t) = 2A2{Q{t)). It is easy to check that 
the condition |||^(t)||| = \\\Ai{n{t)) + A2(fi(t))||| < 2M is valid for any n{t) E S2. By 
introducing /C[^(t;t') such that }C'^{t;t') = W{t)}Cn{t;t')U{t') at t > t' and /C[^(t;t') = / at 
t < t', we can rewrite Eq. (46) as 

^/C^(t;t') = -UHt)[A,{n{t)) + A2{nmU{t)IC'^{t;t') (48) 

when t > t' . The above equation is of the form of Eq. (35), and does not suffer the 
unbounded problem when the Hamiltonian is not bounded. For any two operators ^i{t) 
and VL2{t) G ^2, we have 

sup||/Cn,(t;tO-/Co,(t;t')ll=sup||/C^,(t;t')-^n,(i;OII<2|||fii(t)-fi2(t)||| (49) 
t,t' t,t' 

because sup,>,, ||/C^^(t; t') - /C^,(t;t')ll < sup,>,, || /; rfrf/t(t")[A(fii(t")) + 

A2{n,{t"))]U{r)[IC'^^{t;t") - }C'^,it;t")]\\ + sup,>, II j;,dt"Wit")[A^{n^it")) + A2(fii(t")) - 
Ai(fi2(t")) - A2i^2it"))]U{t")IC'^^{t;t")\\. Based on proposition 3.3, we can prove that 
F(n(t)) G ^2 for any fi(i) G ^2. 

By the following lemma, actually F is a contraction 3^ on ^2 when the product M|t/— to| 
is small enough. 
Lemma 3.5 We can find a positive number c such that 

|||F(fii(t))-F(fi2(t))|||<cM|t^-to|x|||fii(t)-f]2(t)||| (50) 

for any two time- dependent operators Qi{t) and 1^2 (t) G ^2. 

proof Let Ai(t) and A2(t) be two time-dependent operators satisfying 

^A,(t) = z[A,{t),H{t)] - {A,(t), A(^(t))} + {/ - A,(t), A2(^(t))} at t > to (51) 

and Aj(to) = p(^o); where j = 1 and 2. Then Aj(t) = F(rtj(t)), so 

A,(t) = ICn,{t;to)p{to)ICl^{t,to) + 2 f dt'ICn,{t;t')A2{n,{t'))ICl^{t;t') (52) 

from definition 3.4. Here /Cq (t; t') is introduced by Eq. (46). From the above 
equation, we can complete the proof by using Eqs. (29), (44) and (49) because 

\\YiY2Y3-ZiZ2Zs\\ < iiFiiix 11^211 x||y:3-Z3||+||ri||x||r2-^2||x||Z3||+||ri-Zi||x||Z2||x||Z3|| 
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for any six operators Yi, Y2, Y3, Zi, Z2, and Z3. QED 

From the above discussions, the function F is a contraction on 1S2 if we choose tf such 
that M\tf — to I is small enough. Based on the fixed point theory [32], there exists a unique 
pit) E S2 satisfying p{t) = F{p{t)) and we can complete the proof. 

IV. DISCUSSION 

In the last section, a contraction F is defined on S2 such that the unique solution to Eq. 
(16) satisfies pit) = F{p{t)). To obtain such a solution more explicitly, we can construct a 
sequence {pj{t)} by 

Piit)=p{to) (53) 

as t G [to)^/] and 

-p,+i(t) = 2[p,+i(t),i/(t)] - {p,+i(t), Ai(p,(t)))} + {I- p,+i(t), A2(p,(t))} (54) 

= z[p,+i(t),i7(t)] - {p,+i(t),Ai(p,(t)) + A2(p,(t))} + 2A2(p,(t)) 

with pj_|_i(to) = p(to) for each positive integer j. Then Pj+iit) = Fipjit)), and the sequence 
{Pjit)} ^ '^2 converges to the solution we want. 

In definition 3.1, different continuous classes are introduced for the proof. When the 
unitary operator is generated by a bounded self-adjoint operator Hit) G S[, Si and S[ 
are of the same class and it is not necessary to clarify them. However, usually Hit) is 
unbounded, and the time-derivative %pit) should be defined carefully in all the quantum 
master equations mentioned in section II. The meaning of the time-derivative on p(t), in fact, 
becomes more clear if we trasform p(t) to the corresponding matrix p'(t) = f/^(t)p(t)f/(t) in 
the Heisenberg picture. In lemma 3.2, proposition 3.3, and definition 3.4, the density matrix 
pit) G Si and thus p'(t) G iS(, on which the time-derivative can be defined by Eqs. (30) and 

"'■ n 

It is shown in Ref. [8| that Eq. (16) can be reduced to Eq. (6) when p(t) is of the low- 
density distribution such that / — p(t) ~ /. By rewriting Eq. (16) with respect to p®(t), 
we have 

|:p(^)(t) = ^[p'^\t),Hit)] -\Y, ^n'nil - {n\p^\t)\n)){p'^\t), \n'){n'\] (55) 



{n'n) 
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+1 Yl ^n'n{nV'\t)\n'){I - p^\t), \n){n\}. 



(n'n) 

When p^^\t) is of the low-density distribution such that / — p^^\t) ~ /, we can reduce the 
above equation as 

|p(^)(t) = ^[p^^\t),H{t)] - 1 J] c.„,4p®(t), \n'){n'\} + J] uJr,,,\n){n'\p^\t)\n'){n\.{56) 

(n'n) [n'n) 

By comparing the above equation to Eq. (6), we can see that Eq. (16) is also reduced to 
the Markoff master equation of Lindblad form in the low-density limit with respect to holes. 
In Ref. 81], Eq. (16) is derived from Eq. (18) after considering the conservation of number 
of particles in each transition. If we set Ap{t) = 0, we can reduce Eq. (18) to 

d 



dt 



p{t)=t[p{t),H{t)] + {p{t),A,{t)}. 



On the other hand, by rewriting Eq. (18) with respect to p^P\t), we have 

d 



,(p) 



dt 



p^^\t)=t[p^^^it),Hit)] + {p^^^it),A^it)} 



mi 



(57) 



(58) 



if Ap{t) = 0. The above two equations, in fact, can be obtained by considering the non- 
hermitian Hamiltonians with Ap(t) and Ap{t) as imaginary parts for particles and holes, 
respectively. [8] The nonhermitian Hamiltonians have been considered in the quasiparticle 
theory, in which Ap{t) {Ap{t)) can correspond to the lifetime of excited particles (holes) 
above (below) the Fermi energy, [sl, |33|, |3J, |35| The loss and gain factors Lf{p{t)) and 
Gf{p{t)) can be obtained from the last two terms of Eq. (18) by setting Ap{t) = —Ai{p{t)) 
and Ap{t) = —A2{p{t)), so both the lifetimes of particles and holes are incorporated in Eq. 
(16). [8] The master equation, however, may be used to model the time-evolution of the 
nonequilibrium system with no well-defined (quasi-) Fermi energy. Hence the particles and 
holes are the occupied and empty parts of any orbitals in the master equation while they 
describe the filled and empty orbitals above and below Fermi energy in the conventional 
quasiparticle theory. 
The following equation 

|p(t) = z[p(t), H{t)] - \ 5^{p(t), Vi[I - p(t)]V/} + \Y.^I- Pit), V^pim} (59) 

I I 



can be obtained from Eq. (18) by setting Ap 
— \ ^; Vj p{t)Vi. Here {V^} is a set of operators. 



(t) = -lEiHl - Pit)]Vl B.ndAp{t) = 
In the low-density limit with respect to 
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particles (or holes), we can reduce the above equation to the master equation of Lindblad 
form, of which the relaxation term is given by Eq. (7). Eq. (18), in fact, is a general 
equation for fermions and can be used to model different Fermi systems. For an example, 
the density matrix of the form 



P,{t) 



Ps{t) 



Ksit) 

is-pm 



(60) 



has been introduced in the coordinate space for fermionic quasiparticles in the relativistic 
Hartree-Bogoliubov model. 20|, |2l| Here Ps{t) is the one-particle density matrix, K,s{t) is the 



antisymmetric pairing tensor, and Ig is the one-body identity operator. Such a model is a 
BCS pairing model for meson- nucleon couplings. The matrix /^ — p*{t) is the conjugate of 
the matrix Ig — Ps(t) for holes, and Pq{t) is symmetric with respect to particles and holes. To 

preserve the symmetric form of Pq{t) when we model it by Eq. (18), we just need to require 

36| 



SA*(t)S = Ap{t) and S 



/. 
/. 



(61) 



as a constraint on the relaxation term. (In the relativistiv Hartree-Bogoliubov model, the 
effective Hamiltonian H{t) follows SH*{t)S = —H(t).) In Ref. 36|], a similar constraint on 



Eq. (18) is taken into account to extend the time-dependent Bogoliubov equation [25] for 
the quasiparitcles in the conventional superconductors. Such an constaint is important to 
obtain the corresponding semiclassical equation [27] in the incoherent limit. 

In the relativisitic Hartree-Bogoliubov model, the density matrix Pq{t) represents the 
quasiparticles of which the orbitals are composed of electron and hole parts. The quasipar- 
ticles are fermions, so there are corresponding quasiholes. If quasiparticles and quasiholes 
can couple to form "new" particles just as how electrons and holes do, it seems natural to 
introduce the density matrix 



Pq'it) 



(62) 



Pg{t) Kg{t) 

for the "new" particles. Here Iq and Kg(t) are the identity operator and re-pairing tensor 
for the quasiparticle (quasihole) orbitals, respectively. To construct the master equation 
for Pq'{t), we shall consider an additional constraint similar to Eq. (61). |36| We note that 
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multiple order parameters 26|, |28| can be incorporated after introducing the re-pairing tensor 



Hq{t). In addition, both the particle-particle (hole-hole) and particle-hole pairings are taken 
into account in Pq'{t). 36| Two type pairings have been considered to unify BCS theory and 
antiferromagnetic/ferromagnetic theory. 29|, |37| Actually a density-matrix series {pgj(t)} 
can be constructed by |36| 

with j = 1, 2, 3, ... as positive integers. Then we can re-obtain Eqs. (60) and (62) by setting 
Pq,iif) = Ps(t), Pq,2(t) = p J t), and Pq^zit) = Pq' (t) . The corresponding Hamiltonians form a 
chain of Hamiltonians 



PgJ+ii't) 



(63) 



For noninteracting systems composed of finite identical fermions following the modern 
quantum mechanics, the trace trp{t) of the density matrix equals the number of particles 
and thus < trp(t) < oo. It may be reasonable to assume that H(t) = Hq as t < to, and 
the system is in an stationary state at to- In such a case, we shall set 



p{to) = 2_^Cn\n){n\ with < c„ < 1 for all n 



(64) 



such that the initial density matrix is incoherent with respect to the eigenorbitals of Hq. 
Here we require that trp{tQ) = ^„c„ equals the number of particles. Each c„ is the initial 
occupation number at orbital \n). The validity oitrp(t) = trpitn) is expected under Eq. (16) 



because the conservation of particles is taken into account. [6|, |8|, |10| To prove the invariance 
of trp{t), we note that the solution to Eq. (16) is fixed under the contraction F introduced 
in definition 3.4. For any Q{t) G ^2 with sup^trQ{t) < cxd, it is shown in Appendix D that 



trA{t) = trp{to) + J2 ^nn' I dt'{n'\^{t')\n){l - {n\k{t')\n)) 



(65) 



(n'n) 



LUri.T).' 



dt'{l- {n\n{t')\n)){n'\K{t')\n') 



to 



if A(t) = F{Q(t)) when p(to) following Eq. (64) is of a finite trace. If trp(t) is bounded, we 
can set Q{t) = p(t) = F{p(t)) = A(t) in the above equation to prove trp{t) = trp(to). To 
see that trp{t) is finite, we can ignore the third term of the above equation to obtain the 
inequahty 



< trA(t) < trpito) + ^ u;„,„/ / 

(n'n) ■^*" 



dt\n\n{t') \n') < trp{to) + - sup trfi(t) (66) 

2 t 
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since every operator in S2 is positive and we choose tf < to+l/4M. Based on the above equa- 
tion, we can prove that the trace of each Pj{t) constructed by Eqs. (53) and (54) is smaller 
than 2trp(to) by induction. Therefore, we just need to note < trp{t) < lim ^trpjjt) < 
2trp(to) to complete the proof. 

V. CONCLUSION 

In conclusion, the positivity and Pauli's exclusion principle are both preserved under the 
nonlinear quantum master equation introduced in Refs. [8| and [10| when there exists an 
upper bound for the transition rate. Both the loss and gain factors of the equation induce 
the decoherence, which is important to the positivity and Pauli's exclusion principle. The 
number of particles is conserved if the initial density matrix is of a finite trace. Such an 
equation can be generalized to model BCS-type quasiparticles. On the other hand, it can be 
reduced to a Markoff master equation of Lindblad form in the low-density limit with respect 
to particles or holes. 

Appendix A 

For any two kets \a) and |/5) and a time-dependent operator Q(t) E Sq, both 
Re{a\Q{t')\P) and Im{a\Q{t')\i3) are Riemann integrable on [to,tf] because they are contin- 
uous (real) functions. Here Re{a\Q(t')\P) and Im{a\Q(t')\P) denote the real and imaginary 
partsof (a|Q(t')|/?). Hence J^^t' {a\Q{t')\(3) = J^yt'Re{a\Q{t')\(3)+i J^yt'Im{a\Q{t')\f3) 
is well-defined for any ti and ^2 ^ ['^Oi^/]- In addition, the operator O = J^^ dt'Q{t') is 
well-defined in the sense that {a\0\P) = f^^ dt'{a\Q(t')\P) for any \a) and |/?). The time- 
derivative on each operator Q'{t) G S[ is defined by Eqs. (30) and (31) after introducing 
the intergal on Sq. Such a definition, however, cannot be extended to Si when the unitary 
operator U{t) is generated by an unbounded Hamiltonian H{t). 

The following corollaries can be used to simply some calculations in section III. 
Corollary Al Let Qs{t) = Qi{t)Q2{t) for any two Qi{t) and Q2(t) e Sq. Then 
Qsit) G Sq. {Similarly, if ^'^{t) = Q[{t)Q2{t) for any two fi'iit) and fi'2it) G S[, we have 
Q',{t) G S[). 
proof Since |||Q3(t)||| < |||Qi(i)||| x |||Q2(^)|||, we just need to show that 
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/imi2-»ti II [23(^2) — 23(^1)] |tt) II = for all \a) at any time ti G [to,tf]. By setting 
1/5) = Q2{ti)\a), we have Umt.^tAiQsih) - Q3(ti)]|a)|| < ltmt,^tA\[Qi{t2) - Qi{tiW)\\ + 
Umt,^tA\Qi{h)[Q2{t2) - Q2(ti)]|a)|| = 0. QED 



Corollary A2 Let ^[(t) and ^[{t) be two operators in S[, and assume that we can find 

m(t) and Q,{t) = |J 



two operators Qi{t) and Q2{t) in Sq such that Qi{t) = ^il[{t) and Q2{t) = 4^2 (t) on 



[ti,t2] C [to,tf]. Then %{t) = n[{t)n'^{t) E S[ follows m(t) = Qi(t)f^2W + ^'iW22(t) 

when t e [t,,t2]. That ts, |(^(t)f]^(t)) = {§-^n[{t))n'^{t) + n[{t){§-^n'^{t)). 

proof For any \a), we can define two-parameter kets \ai{t,t')) = Qi{t)Q'2{t')\a) and 
|a2(^,^')) — ^i(^)22(^')l<^) ^s (t, t') G X = |ti,t2] X |ti,t2], which is compact under 
Euclidean metric. Both |Q;i(t,t')) and |Q;2(t,t')) are (uniformly) contiuous such that 
liniAt^oAt'^oWMt + At,t' + At')) - |ai(t,t'))|| = HniAt^oAt'^oWMt + ^t,t' + At')) - 
<^2(^,^'))ll = on the compact domain X. Then for any ket \f3), 

^^{p\mt + At) - nmic^) (Ai) 

= ^(/9|fi'i(t + At)(f]^(t + At) - n',it)) + (^(t + At) - n[it))n',it)]\a) 

1 pt-\-At 1 pt+At 



^^ ,. \Q[{t + At)(y^ dt'Q2{t'))\a) + —m I dt'Q,{t'))Q',{t)\a) 



1 rt+At 



^^ , rft'((/3|a2(t + At,t')) + (/3|ai(t',t))) 

|a2(t,t)) + {(3\ai{t,t)) as At -^ 0. 



Hence 



j^{P\nsit)\a) = m[Q,it)n',it) + n[it)Q2it)]\a) (A2) 

for any \a) and \j3) when t G [ti,t2], and we can complete the proof because we have 
1^3(t) G S[ and Qi(t)^(t) + Q[{t)Q2{t) G Sq from corroUary Al. QED 

(2) 
Corollary A3 Let Sq be the vector space composed of all the mappings Q(t; t') 

from {t,t') G [to,t/] x [to,t/] to bounded (linear) operators such that sup ||Q(t; t')|| < cxo 

and/imAt,At'^o||[Q(^ + At;t' + At')-Q(t;t')]|«)|| =0 for all \a). Then each Q{t]t') e S^ 

is intergrable with respect to t and/or t'. In addition, Q3{t;t') = Qi{t;t')Q2{t;t') G Sq 

if both Qi(t;t') and Q2{t;t') G S^ . 
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Corollary A4 Let Q(t;t') be a two-parameter operator in Sq , which is introduced 
in corollary A3. Assume that 0{t;t') follows 0{t;t') = Oo{t') + f^, dt"Q{t";t') at 
t > t' and 0{t;t') = Oo{t') at t < t' for some Co(t) G Sq. Then the operator 
W{t) = fl dt'0{t;t') e S[, and §-^W{t) = Oo{t) + j'l dt'Q{t,t') e Sq. 

proof It is easy to prove that W{t) is well-defined because 0{ti;t) G Sq for each ti. For 
any two kets \a) and \j3), 

{p\W{t)\a)= f dt'{P\0{t;t')\a) (A3) 

Jto 

dt' {p\Oo{t')\a) + f dt'dt"{p\Q{t";t')\a) 

to Jto<t'<t"<t 

= {f3\{ f dt'poit') + f dt"Q{t'-t")W\a). 

Jto Jto 

Then we have fvr(t) = CoW + j'l df Q{t; t') G Sq because \\ Jldt'Q{t;t')\\ < 
\tf - to\ X sup||Q(t;t')|| and HmAt-.oWHJl'^'^^ dt'Q{t + ^^' ^') ~ Itldt'Q{t;t')]\a)\\ = 
for any \a). Since ||W^(t2) - W^(ti)|| < Its -ti| x ||||iy(t)|||, we have W{t) G S[. QED 



Corollary A5 Let Q'(t) G S[ and Q{t) G Sq be two time-dependent operators 
following -^ri'(t) = Q(t) on the interval [ti,t2] C [to,tf]. If Q^{t) G Sq, the equat 
|fi't(t) = gt(t) is valid on [ti,t2]. 



ion 



at 



Appendix B 

To prove the existence and uiqueness of K'{t;t') satisfying Eq. (35), we can construct a 
family of mappings {Qt'} from S[ to S[ itself such that 

ft 



A(t) = / - / dt"Au{t")n'{t") (Bl) 

Jt' 



at t > t' and A'(t) = J at t < t' iff A'(t) = Qt'i^'it)). Such a family of mappings are 
parametrized by t' G [to,tf]. Each Qt' is well-defined because the integrand —Au{t)^'{t) G 
Sq. Since 

\\\g,{Q[{t))-gemtm\<l\\mt)-^m\\, im 
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{Qt'} is a family of contractions on S[. From the fixed point theory, there exists a unique 
K'{t;t') G S[ such that K'{t;t') = Qt'{K'{t;t')) at each t'. In addition, we can introduce 
K{t]t') just as in lemma 3.2. If the Hamiltonian H{t) is bounded, |||^-ft'(t; t')||| is bounded 
and Eq. (36) can be defined by Eqs. (30) and (31). Although Eq. (36) suffers unbounded 
problem when H{t) is not bounded, we can define K{t\ t') by Eq. (35) after tranforming 
K{t-t') to K'{t;t'). 

Direct estimations yield 

\\\gAn'm\\< i + ^-mmi (b3) 

iwg^n'it)) - Qt'.i^'mw < 2M\t',-t[\ X mmi (b4) 

From Eq. (B3), we can prove that 

sup||ir'(t;t')|| <2. (B5) 

t,t' 

(Actually sup^^/ ||fC'(t;t')|| = 1 because Auit) is positive.) In addition, we have 

sup \\K'{t;t'^)-K'{t;t[)\\ < 8M\t'^-t[\ (B6) 

t 

from Eqs. (B2) and (B4) because K'{t;t') is fixed by Qt'- Eq. (B6) indicates that K'{t;t') 
can be taken as as a mapping from t G [toj^/] to S[ C So, so K'{t;t') is integrable with 
respect to t and/or t'. 

Since K'{t;t') G S[ at each t' G [to,t/], we have /imAi,At'^o||[^'(i + At; t' + At') - 
K'it- 1')] II < /imAi,At'-o II [K'(t + At; t' + At') -K'(t + At; t')] || + /imAi,At'-o II [K'{t + At; t') - 
K'{t\ t')] II = at any (t, t') in the definition domain from Eq. (B6). Together with Eq. (B5), 
therefore, we can prove that K'(t;t') and K''^{t]t') G iSq , which is introduced in corollary 
A3. It is easy to see that K{t] t') and K\t] t') are also in Sq from the definition of K{t] t'). 
From corollary A3, we can prove that the integrals in Eqs. (37) and (38) are well-defined 
because Bu{t') and B{t') can be taken as operators in iSq . 

Appendix C 

In this Appendix, we assume that Vl{t) G ^2. Since AiiVlit)) and A2{yL{t)) are both diago- 
nalized by the eigenorbitals {|n)} of i/o) we have |||Ai(r2(t))||| = sup„ J(n|Ai(r2(t))|n)| < M 
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and |||A2(r2(t))||| = sup„( |(n 1^2(^2 (t))|n)| < M after some calculations. To show that 
Ai{yt{t)) and A2(f2(t)) G Sq, therefore, we just need to prove that 

Umt,^tAMn{t2))\a) - A,{Q{h))\a)\\ = (CI) 

for all \a) as j = 1 and 2. It is easy to see the validity of the above equation if \a) is in the 
set composed of all the finite linear combinations of eigenorbitals of Hq. Because such a set 
is dense and both |||Ai(t)||| and |||A2(t)||| are bounded by M, we can prove that Eq. (CI) 
holds true for any \a). 

Appendix D 

In this Appendix, we assume that p(to) satisfies Eq. (64). To prove Eq. (65), 
in which A(t) = F{Vt{t))^ it is important to estimate trKLQ{t^ta)\n){n\K/^{t-,ta) = 
{n\ICl^{t;ta)ICn{t;ta)\n) for each normalized eigenorbital \n) of H{t) because both p(to) and 
A2{fi{t)) in Eq. (45) are diagonalized by {|^)}. Here we assume that t^ <ta < t. Because 
(n|4(t;t„)/Cf,(t;t„)|n) = (n|f/t(t„)/C;^(t; t„)/C^(t; 4)f/(4)|n), |(n|4(t; t„)/Cn(t; tj|n) is 
well-defined and equals -2{n\K\^{t]ta)[Ai{VL{t)) + A2{VL{t))]lCn{t-,ta)\n). So 

tr/C^(t;4)|n)(n|4(t;tJ = I - J2 [ dt'\{n'\}Cnit';Q\n)\'' (Dl) 

x[un"n'{l - {n"\n{t')\n"))+Un>n"{n"\^{t')\n")]. 

Let a (time-independent) operator C = '^nCn\n){n\ with {C„} as a set of positive real 
numbers such that J^n^n < oo. Such an operator is diagonalized by the eigenorbitals of 
Hq. We can extend the above equation as 

trlCnit;ta)CICl,it;ta) = trC- J^ [ dt'{n'\ICn{t';ta)CICUt';ta)\n) (D2) 

x[ujn"n'{l - {n"\^{t')\n")) + uJn'n"{n"W)\n")]. 

Let {Cj} be a set of finite (time-independent) operators parametrized by j such that 
Cj = Xln^" l^)(^l- Here {Cn } is a countable set of positive real numbers such that 
X]„C'„ < oo for each j. Let 

r(t) = Y, trJCnit; t«)C,4(t; t,). (D3) 



t,<t 
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We can prove that 

tj<t {n'n)"'*0 



trT{t) = ^trCj -^ f dt'{n\T{t')\n)[uJn>nil - {n'\n{t')\n')) (D4) 

t .^i r..i.^\ J to 



+uJnr.'{n'\nit')\n')] 

from Eq. (D2). Here each time tj > to- Actually Eq. (45) yields a form similar to that 
provided by Eq. (D3) because F{il{t)) equals /Cn(t; to)p(to)^Q(^; ^o) plus an integral term, 
which can be taken as a limit of a summation over }CQ{t;t')A2{i^{t'))}Cl^(t;t'). Because 
sup^trA2(f2(t)) < M X sup^trQ{t) for any Q{t) G ^2, we can generalize Eq. (D4) to prove 
Eq. (65) when both trp(tQ) and sup^trQ{t) are finite. 
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